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I. INTRODUCTION 


The idea that both atmospheric and subsurface environ- 
mental conditions can have a Significant =ffect (enhancement 
er degradation) Oh modern weapon and sensor systems has now 
become a reality to many military leaders. This has led to 
a greater demand for more reliable and accurate atmospheric 
and oceanic prediction schemes. Prediction of the existence 
and evolution of electromagnetic ducts are important in pre- 
mening Over the horizon surfaca to surfase detection. Pre- 
Sa Ce 


= 
— 


dictions of the inversion heijght is important wher 


(D 


Known that optical propagation is degraded due «=o zurbulence 
and extinction. Ocean mixed layer predictions are importan+ 
mee cdete=mining convergence zons distances, optimum source 
Memes Gep th, and cutecrr frequency for acoustic surface duct 
propagetion. 


The abiiity to rrovide accurate forecasts requires 


a) 


Sempiets undsrstanding of the dynamical orocesses occurring 
*-n each of the atmospheric and oceanic boundary layers, and 
as important, the processes 2s5sociazed with coupled iccal 


charges in the wo layers. 


Venyebitceie is Known about the interactive nature of the 
edjacent turbulent boundary layers. Pee ogeectave Gf chs 


i 





two already 


thesis is to formulate a model which couples 
proven models for the respective adjacent 

This coupling is achieved by 

stress, mass 


existing and 
layers. 


well-mixed boundary 
surface boundary conditions for 


used to examine the pro- 


matching the 

flux and heat flux (Figure 1). 
This coupled model may then he 

and atmosoheric mixed layers 

hours and 50-100 ka, 


cesses of coupled local oceanic 
al time and spva*ial scales of 12-18 
respectively. Sede ctGadkht Changes in both mixed layers can 
these scales changes requires 
oceanic 


oa 


Explaining these 
Sivexisi ng Acmesphsric and 


mecur at 
coupled evaluations 
well-mixed layers at a point. 
< the present time, we only have qualitative evidence 
“hat changes in each layer ars correlated. The figcrous 
Seem csGeioe ONS Sf the coapiing will require 
CESOe cS 2nd, perhaps, improved 


eee Ciricatior 
considerable interpretive 
formulation of the separate moiels 


in 













Begure 1. 









E/M PROPAGATION 
MODELS 


ATMOSPHERIC initial 
ANALYSIS aa ae 
conditions 
as, ATMOSPHERIC BOUN- 
DARY LAYER MODEL 
SATELLITE, wind heat, moisture 
A/G, SHIP stress fluxes, 
OBSERVATIONS 
ae 


OCEAN THERMAL 
STRUCTURE ANALYSIS (3) 


radiation 
= Je 
OCEAN MIXED LAYER 
; MODEL 
initial | ae ee 


condi- 
tions 







OCEAN GCM (1) 





SOUND VELOCITY 
SIRECTURE 


ACOUSTIC PRCPACATION 
MODELS (4) 





Coupled Environmen=2l Model System. (1) Embedded 
mixed layer in OCGM eae Seal, 1981) OF 
TOPS (2) NCGAPS/NORAZS or 1-0 nodsl with 
prescri ed subsidence (3) EOTS or other (4) At 
BuegGe Or TC APS/SINIS,- ¢ct. 





II. BACKGROUND 


A. GENERAL BOUNCARY LAYER FEATURES 

The Marine Atmospheric Boundary Lay2r (MABL) extends 
from the surface *hrough the capping inversion. The Oceanic 
Boundary Layer (OBL) extends from the surface to the sea- 
sonal thermocline. The air-s3a interface is bordered by 
oceanic and atmospheric turbulent mixed layers which effec- 


ely insulate the quasi-geastrophic  %ceanic and atmos- 


ct 
< 


2 


pheric flows below the thermocline and abov2 the inversion, 
respectively. The primary source of hs turbulence is the 


velocity (current) and buoyancy (density) gradients created 


(Dp 


Byeanhteraction a+ the interfaces. Mhewweege  “vertscal mixing 
yields nearly homogenscous (weli-mixed) wind, temverature and 


Mmimerdt.y profiles above the surface. In the ocean, the *ern 


t(D 


Mmmm 14 Yer rerers #£c tne region cf vertical uniformity in 
mean velocity and density. ALt howe@h density ss a ftunctzen 


Semooch =empexzature and salinity, it is usually d 


(Dp 
Tw 
Mm 
Qu 
D 
ee 
ct 


upon temperéeturs except when large near-surface 


i) 
ey) 
f+ 
|- 
: 
| 
ct 
ht 


Sie ges occur due to heavy precipitation or ice melt in high 


{+ 


la-+ 
aaa 


tudes. At the bottom of the ocean nixed layer and at 
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the top of the atmosrhere mixed layer are thin “<ransition 
tegicns (thermocline cr inversion). The mixed layers inter- 
act with the upper and lower layers at these regions by 
means of turbulent forced entrainment. Typical model pro- 
files deépicting the boundary layer regions for both <«he 
atmosphere and ocean are shown in Figure 2. 

The weli-mixed nature of both the turbulent MABL and OBL 
mip lies features for the vertical distribution of mean val- 
ues and vertical fluxes of wind, temperature and humidity in 


the atmospnere and density and meméentum in the ocean. One 


t+- 


feature is that the properties which ar2 conserved during 
mixing can be treated as being constant with height (depth) 
in the boundery layer. A second featurs is that vertical 
tluxes of the well-mixed parameters vary iinearly with 
mei gnt. These features enable predictions to be based 
solely on the surface and inversion (entrainment PLU Kes we 
advection is negligible. 

Maer ner UMporttant source of energy, in addition to tur- 
bulent kinetic energy, Pieoern Maced Layers: iS radiation. 


Mera ke the MASL, most of the solar radiation does not pene- 


Ih 
4 
= 
rs 
3 


trate the OBL. Thererore, downward curbulernt heat 


the OBL is aS important as the upward flux. 2hee whe long 
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emitted 


flux is absorbed and 
Sltewce=e 0 this flux 


radiaticn 
~~ 


of the sea 
the oceanic 


(infrared) 
m2 1limeter 
componente 


wave 
fest 
of 


essentially anoth2r 
surface heat flux. The total surface heat flux is comprised 
heat and the long and short wave 
(solar) 


and latent 
and short wave 


Eien Lo ng 
dependen*= upon 


wethin ihe 


component is 


che 


of sensible 
Bcth 


madea- lon. 
Gare ucaliy 


SGC en vare 


solar 
into the 
buoyantly-creaced velocity 


fluxes 


existance of clouds in the atmosphere. 
on 


and OBL, 
effecz 


MABL 
have the most proneunced 
ends 


in both the 
Peuctuations and mixing 
mixed layer characteristics. Garwood (1977) states 
vertical mixing has a nore obvious and 
energy budget, as shown 


buo yancy-irorced 
direct role in the mechanical 
3. In the atmosphere, 2ven under near 
Peevcr is «SOPs sand be quits 
eddies 


Energetic 
Toes he invews: oO 


= 
ae 


Fiqure 
buoyercy forced 
driven 


These buoyancy 
the surface 


eons, 
fron ch 
more nome 


Reum 


large. 
atmospheric boundary layer 
day Sait bean ee 
causes 


and 


Moe-e ~hey entrain warn, 
mixed layer. Tia ne 


e1trainnent 


Laver 


=nto the 
excend above the lifting ccndéensation 


feet TOrm Within the layer. 


ie 
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; PLUS ABSORPTION OF note 
MEAN K.E, 
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HEAT 


Figure 3. for the Ocean 


Tt is the realization that buoyancy driven 


depends on the near surface stability, wheeh “on 


Bererclied by local atmospheric forcing >£ «he 


dayer which makes a coupled pr2diction approach nécessary. 
MeeeeecCUularity in cause and effect ris2s because surface 
Doyency flux induced entr2inmen= not only increases the 


Mee Of the MABL, but it also changes its effect on the 
ocean mixed layer. TAS Bele DE the surface buoyancy fiux 
means changes in the oceanic well-mixed layer could be 
€ither the cause or the resuit of entrainment. Another 
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cause and effect relationship is the forma*ion and dissipa- 
meon Of clouds within the ABL. Clouds can be caused by 
changes of the ocean surface tSmperaturse. iie=ten, clouds 
have a profound effect on the short and Long wave radiation 
budget for the OBL. 

These cause and effect ralationshinos shew that the 
assessment of rélative roles of dynamic processes in both 
layers crequires coupled evaluations using relatively compi¢cx 
Maysical models. A data set collected (Davidson e= al 1983) 
ever a 48 hour period in May 1978 demonstrates *he nature of 
coupled local changes in the atmospheric and océanic bound- 
eimy layers. MObma weer loa macm 1700 5 0500 local ¢imes, 
clouds 2xisted in the MABL and winds increased steadily from 
6 to9 ms! duGpane@wenis “period. Te is suggest=d czhis 
encrease was associated with deepening of the MABL (and 
Smerercore, the entrainment of additional momentum into «he 


fexed laver) which is, in turn, associated with the surface 


iD 
+e | 
PD 
fu 
fv 
WV 
ct 
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temperature. The oceanic boundary layer deseo 
Wind reached @ maximum, tesuiting in the surface <emperature 
m=OUG.ng rrom 14.59 £0 13.09C. 

Theré were other atmospheric tespons= and forcing feae- 


Baees Observed during the entire 48 hour period from 


kee, 





5/19/1200 to 5/21/1200 PST which are balieved tec have been 
associazed with coupling between the two layers. Atmos- 
pheric and oceanic mixed layer observations appear in Figure 
4. 

These additional features are: 

a) The ABL depth increased from 250 to 750 meters over 
the 48 hour period. The changes occured in relatively 
Sore shrtomvars, £rom, 2O7 0000, 25 20/0400, and fron 
21/0000 to 21/0500. Tha ABL depth remained at near 
constant depths during the intervening 12 to 18 hour 
periods. 

b) The surface layer temperature increase from 20/0000 
momces7 i500 2s sindicative of entrainment of overlying 
warm air. This entrainnean* presumably was a factor 
in the wind speed increase over the same period. The 
speed increase, in turn, presumably caused <he wara 
Shallow OBL to be destroyed which dropped «he surface 
temperature after 20/1800. 

€) The most dramatic chang in the MABL occured during 
the period from 21/0000 to 21/0500 which was immedi- 
ately after «he surface temperaturs had dropp=d due 
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Questions which could be posed on the basis of the above 
data are: 

a) Did the dramatic change in the ABL depth at 21/0000 
occur as a result of the ocean forcing the atmosphere? 

b) Does the fact that the mixed layer evolved with rapid 
changes followed by equilibrium periods indicate fesd- 
back between the boundaries? 

c) What would have occurred if the wind had reached its 
Maximum value at a time Other than after suns2t? 


d) Was the time of occurrence of the wind maximum and, 


(Dp 


hence, «he ocean mixing controlled by the coupled 
mixed layer systen? 

This thesis will attempt +o answer some of these ques- 
mons On the basis of the formulated csupled model. As 
stated earlier, all answers can only be obtained through 
considerable interpretive efforts. 


The scope of the overall scientific ffort is illus- 


(D 
\3 


mee ed 1n Figure 5. RicwiOeemmmOcecans CC G=SCr = ptions will be 


made by a bulk oceanic mixed layer model (the Garwood 


Model). The atmospheric descriptions will be made by «he 
Naval Postgraduate School (NPS) Marine Atmospheric Boundary 
Layer model (Davidson et al 1983). Pactweweqey Szgniticant 
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operational descripticns (forecasts) to which the improved 
understandings would apply, appear on the extreme right hand 


side of Figure 5. 


B. MARINE ATMOSPHERIC BOUNDARY LAYER (MABL) MODEL 

The NPS boundary layer mod2l is a zero-order, two layer, 
integrated mixed layer model. The two layers censist of the 
well-mixed, turbulent boundary layer, and the relatively 
non-turbulent free atmosphers. The model is based on 


entrainment energetics formulated by Stage and Businger 
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(1981) and cn radiative transfars describe 
oe (1983). 
The inversion (or transition zone) separates one layer 


mon the other. Tn a zero-orier nodel, -nelinversi on os 


(D 


meoumed tO have zero thickness, hence, 2 discontinuity, or 
Sp OCCULS ir the profiles of the conservative parameters 


MVELELON. 


1+. 


at «he 


The present model requires as input, an initial atmos- 


pheric soundirg, the mean winds at a level within the sur- 
face layer (10-30 meters) Semen e SUStace “sempehature. 
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Beundary layer wind values for the duration of th 
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period can be obtained from larger scale models. 


Scudy, actuai hourly wind spe2sd measurenenzts were input as 
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Pees Of the model initialization. Sea-surface temperature 
remains unchanged from the initial input value. However, in 
the coupled mcedel, the sea-surface temperature will change 
during the forecast period and will be that predicted by the 
OBL modc2i. Well-mixed temperature and humidity are predic- 
+ed soc the surface wird and wind shear at the inversion are 


the only atmospheric variables which have +o be prescribed. 


'y 


The larger scale subsidence, cbtained from the large scale 
modé¢l must also be prescribed for the forecast period. 


Three methods can be used tO compute the large scale 


retical velocity (sultsidence) from single station observa- 


< 
Dp 


} +- 


ons. These three methods are well described by Gleason 


fc 
fro?) and include the kinematic and adiabatic methods end 


q 


"he integration of the moisture budget equation (Q-metheod). 


Gleason's study showed the Q-method displayed he most mnerix 


as a Single-station assessment of subsidence. This method 
Was used to compute subsiderc2 for this study. Th] ‘solas 


menich angle used in the computaticn of =he short wave radi- 


(1) 


meeve G1NX is determined from ‘the latitujie, julian day ane 
meme OL Gay of the initial input. 
The model predicts the time evolution of the inversion 


height, and the mixed layer values of 9 (:emperaturs) and gq 


ZS: 





(humidity), as well as their jumps, at 30-minute intervals. 
The predicticn model is run out to 24 hours. The steps in 
the predicticn computation are shown in FPigure 6, where it 
is noted tha+ procedures are the same for clear and cloudy 
cases except for entrainment computation and estimation of 
Siloud top cooling. 

The model uses the standard integrated rate equations 
(Ternek2s and Dreidonks, 1981) to predict the time evolution 
ct the conservative variables and their jumps at the inver- 


Sion. These equations are: 


h(Dx /Dt) = Hee), -(w'x'), + source (1) 


h (D4x/Dt)= hfx (dh/sdt) - rere) a WX) = Sietiece (2) 


where [ is the lapse rate above the inversion and the source 
ae LS €quai to -(F,- ia oes =a eqwal ©O ZSr0 FOr 
Pere oS che net radiative Elux. The subscripts zero and 
h refer to surface and inversion height values, TSSpec- 
fayely. 

cage and Businger’s (1981) entrainment velocity parame- 
terization is used to close the system of equations and 


determine the time evcluticn of the inversion height. The 


cy 
th 
ck 
=m 
—) 


Seostte aSsumpition states that the dissipation rate o 
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TEMPERATURE, HUMIDITY 
SURFACE TO 3 KM 






PRESURUGE. 
SURFACE WINDS 
= SEA SURFACE TEMPERATURE 
SUBSIDENCE (DIVERGENCE ) 
ADVECTION 


COMPUTE : 
SURFACE FLUXES 
CONDENSATION LEVELS 


Gre"s0 ) NO clouds YES 


COMPUTE : 
ENTRAINMENT 














COMPUTE : 
ENTRAINMENT 
CLOUD TOP COOLING 
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moguse 6. Schematic cf Input, Prescription and Computing 
Seees on ap. Predict ion. 
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ere ouient kinetic energy (TKE) is a fraction {1-A) of th 


WM 


preduction rate. WAeeersce che Strat nmen= coefkficzent, taken 
er 0. 2. 

Mathematical modeling of radiative flux transfer is an 
extremely complicated but still developed subject. Unfortu- 


nately, even the simple radiative transfer models are stili 


4 


extremely complex in comparison to most other physical 
parameterizations used inthe mixed layer model. Uncer- 
tainty in background aerosols, atmospheric absorbing gases 
(water vapor, carbon dioxide, »zone) and cloud droplet size 
Pectta are sources of error in radiation calculations. 
Long and short wave radiation fluxes are computed sepa- 
masaely. 

The long wave radiaticn fiux calculation was modified <o 


permit non-black stratus clouds by introducing th CLOuE 


(D 


emissivity, € , which is a function of the total cloud ligq- 
mea content, Ww. Br mcewerewcl oud Liquid water orofiles are 


approximately linear with height (Davidson e+ al, 1982), 


cloud iiquid water content ani clioud emissivity are aqiven 
Oy: 
We 0.594 (h ~ 2) di, (3) 
Eo = 1 - exp (-aW) (4) 
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where Pa _omene density Of a@ar ( 1. 25X 10°% gm i Se 1 
Meeght Of the mixed Jayer (cloud top), Ze is the lifting 
eondensation level (cloud bottom), a=0.158 (Slingo, 1982), 
and q,, Bometre iigu'd Water Conten= at th2 cloud «op. 

ime tong Wave Cloud top nec rediation flux, L,,, can be 
calculated from the cloud top temperatura, The and the 
peeecetivc tacd-ative sky temperature, To,.,, using the Stefen- 


Peeeeztann Law. At the cloud bottom, this fiux, L,., is cal- 
nc 


culated in “he same Manner using the sea surface 
Semperature, T,, and cioud bottom temperature, T,. These 


fluxes are given by: 


eee an ~ Bais) (5) 


Lac nS = tT. ) (6) 


Sse GC -S Stefan's ccnstant (4.61 X oe ) ana. ,  ‘Cuxeid 


Piesclivi-y, is obtained from equation 4. The net long wave 
peux a= che surface, Bera , Decomes: 

= eS 4 =4 4 

Pionges C-s0 Cee 7 (1 €,) Tsky ) (7) 
Phere Tis <he averaqe temperature of the cloud. NOs he 
@moud free case, milemiCeurlux= Ss eleecalculazred at Z=h and 


Z=0 uSing the standard method (Fleagle and Businger, 1978) 
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Semeencegrating the £lux emissivity profil 


Meee citxX at the surface for <he claar sky case is: 


29 





= F = er (8) 


Fiong u 


where F, is the upward radiative flux and Fy ss the downward 
medi:ative flux. 


Surface flux parameterization (bulk asrodynamical) for- 


muilae 
1 *= ona (2)) 
T= Cy (89-9) (9b) 
1/2 
Ge Cy (dos q) (9c) 


are used to determine the surface fluxes »f momentum, sensi- 


bie heat and latent heat. Thesa fluxes are given by: 


2 
u'wi= ux (momen tum) (104 
Tiw'= uu *t* (sensible heat) (10b) 
q'w'= u*q* (laten= heat) (10c) 


Meee Cy and Cy are the ten meter stability dependent drag 


S@eerlclent, 9 iS the potencizl temperaturs, andqiis the 
Peecitic humidity. The subscript zero jienotes the surface 


value. 
Mie shoh= Wave radiative flux is calculesed using the 
Gdelea-fdingtcr method (Joseph, 1976). This was added to the 


Meme l => accoun= for heating of the nixed layer by solar 





Bama Ome) SnOr. wWaVe EXtinetion is dominated by scattering 
(as opposed «c absorption). tiioeseaeces sa radiation forms 
a second short wave radiative component usually referred +o 
meee datfuse® solar radiation. Scattering in the mixed layer 
is due to atmospheric particles, cloud droplets in the 


cloudy case and sea-sait aerosols in the cloud-fre¢ case. 


| 4 


Ffuss radi 


{- 
a 
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This MABL model computes both the direct ard d 
tive components *o determine th total short wave radiation 
flux at the surface. An explanation of the delta-Edington 
methcd and all of the parameters, atmospheric factors, and 
equations involved is a very complex subject and is beyond 
mae, scope of =hiS paper. An excellent review of this sub- 
ject has beer published by FPairall (1981). An 2NpOs-an= 
short wave radiative parameter that is prescribed in the ABL 


model is the fraction of reflected short wave radiation, Ag, 


tty 


neon th ace. Tet iS taken as 0.1. 


(D 


pee Ue 


es OCEANIC BCUNDARY LAYER (OBL) MODEL 

Gartwocd (1977) developed an ocean mixed layer model 
USing the Navier-Stckes equation cf motion with the geostro- 
phic component eliminated, the continuity equation in incon- 
pressible water, the heat Squation fron the first law of 
Meee modynamics, the conservation of salt equation, an ac 


linearized equation of state. 
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PeccmecttemmeOr tle Eat] Of dSepening (sc tetreaz) of the 
mixed layer is dependent upon an understanding of the dynan- 
mes Of the entrainmen*= process. The turbulence of the cover- 


lying mixed layer provides the energy netded to destabilize 


and erode the nderlying stabl2= water mass (Garwood 1977). 


Therefore, «he turbulent kinetic energy budget is the basis 
for the entrainment. This system is closed using 2 mean- 


-? 


turbulent-fieid modeling of th> vertically integrated equa- 


a 
\ 


O@n= for the individual TKE components, plus the inclusion 


@f the bulk buoyancy and momentum eguaticois. 


tal 


Seperate vertical and horizontal equations for IKE are 
used to better handle the mixing process. Buoyancy produced 
energy is somewnat nore efficient than shear production as a 
Bemmce SL energy for vertical mixing because of its uniqus 


Gam ecOnponcit or he =e bulen= velocity. 


pv 


> 


& = oe ae ot 
Serccs Or the vert 


l 
1 


The bucyancy equation iS generated from the heat and salt 


equations together with an equation of st2+2: 
P= Po fl-a@ - 6) + B(S- sy (11) 
mam puoyancy is given by: 


DE gf fo - P)/ Po (12) 


BV 





where @ is temperature, S is salinity, 9 Veo oy,) de Grave ty, 
and the constants « and 6 are t=he expansion coefficients for 
heat and salt, respectively. The *ilde represents the total 
instantaneous value and the subscript zeros denotes an arbi- 


: Using b instead 


th 


@eery, but representative, constant valu 
cof 9 (temperature) EO@MmeneCeP oS fo na con of bucyancy allows 
this model to bé applicabie in situations where evaporation 


Meee cipitation contribut> significantly to the surface 


O 


Biwoyancy flux. Temperature is in most cases considered ¢ 
Bemtee don ne=ing factor in the density vorofile. The rela- 


“ive importance of the salinity effect on the short tern 


emt on Of =he density profile is generally non signiri- 
eamt (Milier, 1976). 

Temperature and salinity profiles, and for extended 
mee@ecasts, the wind driven nhozizgecntal current profiles, are 
merge cd as input for model initializatis2. These are used 
Semcoempute the mixed layer depth. The Sazrwocd OBL modéil 
defines the mixed layer depth, h, as *he shallowest depth 
@emwnich the cbhserved density value, a; , SS Oc ors Unt Ss 


greater than the observed surface density value. 


t 


Ocean envircnmental parameters that have +o be ore- 


ecribed include; theme aee=o OF “Short Wave radiation 
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absorbed in the top one meter of *he ocean, the radiation 
extinction coefficient for datermination of *he downward 
mui bulen: heat flux, and the critical Richardson number for 
a stability adjustment at the bottom of the mixed layer. 

The parameters required for surface boundary condition 


computations include wind speed and direction, cloud cover, 


™ 
rh 


sea surface temperature, air ztamperature (dry bulb), dew- 


Beer. =~smperature, incident solar radiation, and the rate of 


eyroperetion (E) and precipitation (P). 
Pee cturbuien= tluxes ci laten= heat, 0,, and sensible 


[Pete o., 42 €Stimated using the bulk 23rodynam-c formu- 


th 
ij 
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mmemthne Net back radiaticn is estinatad 


eera-20n (Husby and Seckel, 1975): 
-7 4 
eee 8e 4X) 10 (273. 164+T,) (.39 - .5E 


Mteeest. -_S the Saturated vapor voressure of the marine aiz 


Meese cotrec=s fer salt defects), ro 
Sieear based on the dew-point temperaturs, TT, air tempera- 


mec, Tf. Sea surface temperature, Stee cone taacesere | 


meeud Cover. The upward hea= flux, Q,, is =hen given as: 





Cjamest Clb. (14) 
Wame-olar fbedtats.on, O., is given by; 
b 3 

Opa oa oes) 0, (15) 


mere 0, is =he clear sky radiation given by (Séckel and 


Be2uday, 1973): 
G@,= Apt 4,ces Q + B,sin Q + A,cos 29 + B,sin 2k (16) 


The constants a and b are adopted from Tabata (1964) and the 
emmec cloud cover ccrrection from Lasvastu (1960). La 
variable is the mid-day elevation angle of *he sun. ine 
Peeeeetc ents (A,,A,, €Ct.) were calculated by harmonic rep- 
resentation of the values predicted in *he Smithsonian 


Sc, 1958) with 


i-! 
| be 


Meteorologicéel Tables ( 


Q= (277/365) (*- 21) (ay 


where * is tne julian day of the vear. 


Roc all of the incomin Short wave radiation venetrates 
me ocean mixed layer. Approximateiy half (for oven océan) 


Memmeapsorbed within the first meter. The amcunt abscrbed 


() 


MeeeeS from =zegion *o région depending on the amount of 
absorbing particles such as phytoplankton and yellow sub- 


stance. More radiation will b3 absorbed in coastal regions 
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than in the cpen ocean regions. Mins S8Onc2On Cf absorbed 
short wave radiation is therefore consid2red to he part of 
the upward heat flux. The ramainder of the short wave radi- 
ation does penetrate the mixed layer and is attenuated in an 
exponential fashion depending on water turbidity. Tf the 
mea@ctson cf Solar radiation absorbed in the first meter, Rf, 


*s given, ther the net heat at the surface is given by: 


Cet eo came ee (18) 


From these, the surface fluxes of buoyancy (heat and 
salt) and momentum can then be computed. The turbulent ten- 


Pe=atUrs, séliznity, velocity and buoyansy fluxes are given 


Dy : 
(Thw') = -Qnet 7PCp (19a) 
(S'w') = (P=E) So (19b) 
{u'wi) = ux’ (19c) 
(btw!) = ofa (Ttw') - G(S'w')| (194) 
mespectively, where u*, friction velocity in ait, is calcu- 
lated by: 
c= ese 20 
Ss We 10 ( ) 
12 
us= (T./P, ) (2%) 
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The subscript zero refers to ths 


-2 
surface stress (dynes cm §). Ee 


buoyancy resuits when Qnet<O and EDP. 


ods, negative buoyancy flux resuits 
the solar heating at the surface, 
Beriods, a positive 


bined effect of net leng wave 


bulent fiuxes cf heat and moisture. 


New ocear temperature, salinity 


ProttiesS are predicted at one hour 


then used to predict the new 


miene QoEediction computation are shown 


D. FORMULATED COUPLELT BOUNDARY LAYER 


The physical processes described 
Sections were not altered during 
Moet: cat~ijons were involved on 


between the nodels. Also nes2essary 
units. 
use in atmosptreric computations 
temperature. 


atmospneric paramet 
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Se 


while 


buoyancy flux is produced 


ama wind d=iver 


mixed layer i 


— - cell ome tg 
eensi2cs ing 


The only parameter supplied fr 


HOWEVGE, ches Passa lection is 


surface value and @, is the 
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Parameters in the atmosphere 


wind speed and direction, 


julian day, and time of day. 


(direct and diffuse), iong wav2 ra 


jatent heat fiuxes, and the surfa 


puted by the atmospheric section 


-nput <5 the ocean model. 
The first problem associatid 


me differen: prediction ‘time step 


~ 


The atmospheric medel uses 2 


jv 


whereas, the ccean model uses on 


Assuming neo Significant sea-surfac 


take place over the short one-h 


el's prediccion scheme was inse 


model 2s a subroutine. Jel mala! en’) 


called every other 30-minute tine 


eae = Ke eu 


eens. n <¢s 


The second prebiem was wit 


does not require wind directica as 


“he magnitude of the wind speed. 


meguares Wind direction to compute 


Bo 


dew-point temp2rature, 


moijel's 


wind 


Mites lL ization -includ=:: 


eee aS, 


Short wave radiative flux 


diative flux, sensible and 


ce momentum flux were con- 


Se en] Coupled model and 


with the coupling involved 


S 
30-minuzs 


e-hour 


€ temperature change would 


+he ocean nod- 


rted che 


atmospher> 


e s 
es as 


— <—<_ wf © 


= eid 2 scheme only 
Step. 
meg 


Wiktesa Le zat LOR. 
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These components are used to determine the 

Peo 
(Gut, COMpPoOrents used an the 
A new Subrodutine was added 


ee che winds. 
Gcean turbulent velocit 
in 


momentum pudget computaticns. 
scheme to 


@me initialization determine horizontal wind 
components frem wind speed and direc These conm- 


+ermined by: 


ponents are deter 
Boy = 2 (S22 9) 445 (22a) 
Uioy = ~(COS 9) Y4Q (22b) 
WhElTE U,p 1S the specified wind speed at 10 metezs and the 
= A 
COMPChENTS ATS Winx 26S Uigy ¢ 


and north-south 
Winds are 


From which the 


east-west 
respectively. Seis ecie d2tece.on Fr 
Peewiang relative to true north. Use of squatiorns 20 and 21 
and the following relaticnrship: 
2 
Us= Pal*® = Pyiry (23) 
Er Eoeevelocity flux (nt) 20 BS Cetes> 
gS 2 bie ee 


the ocean 
uming the density of Seawater, py, 


Ing the 


ailows 


mined (ass 
mated by 1) Dye OES tees Mods = 
atmospheric parameters and «hen eauating 

wi tress value. Pas fons 


vsiocity flux *o the computed wind s 
SMe OLS 


Volse ray fluxes Ps,, and 


ZOntal components of the 


2 _ 2 
aye, ate ther easily computed by: 
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2 2 2 2 
eM ee Qltox se oy ? (24) 
2 


Wty y2= Cy tyoy (25b) 


The atmospheric derived values are computed using the MKS 
iio system and had to be cConveartei to the CGS system prior 


s done at zhe 


momany Oceanic computations. fhis conversion 
Beert of each ocean prediction run. 

As discussed in section B, the short wave radiative fiux 
is calculated using the delta-Edington method. The long wave 
radiative, sensible heat and latent heat fluxes are computed 
Meeeng equaticns t0a, 0b and 410c, respectively. A unit con- 
Zeasion Bron Watts m <*to calcm s 1S S-dilesed sor these 
values prior to use by ocean model. This conversion is 
mgeen done a~ the start of the Ocean prediction run. A flow 
diaqram of tne steps in che coupled model's prediction cen- 


Biee~20n is shown in Figure 8. 
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mecqure 8. 








INITIAL PROFILES: 
ATM- TEMPERATURE 
HUMIDITY 


OCN- TEMPERATURE 
SALINITY 








Call Ocean Model 
Compute New: 

Sea Surface Temperature 
Mixed Layer Depth 








PRESCRIBE: 
SURFACE. WINDS 
SEA SURFACE TEMPERATURE 
SUBSIDENCE (DIVERGENCE) 
ADVECTION 






COMPUTE : 
SURFACE FLUXES 
CONDENSATION LEVELS 


> NO clouds YES 


COMPUTE : 
ENTRAINMENT 
CLOUD TOP COOLI! 









COMPUTE: 
ENTRAINMENT 







COMPUTE NEW: 
MIXED LAYER DEPTH 
WELL MIXED TEMPERATURE 
AND HUMIDITY 
JUMP STRENGTHS 


COMPUTE NEW: 
EM DUCTS 
OPTICAL TURBULENCE 
AND EXTINCTION 
DISPERSION 





schematic of Atmospheric. and Oceanic Coupled 
Model. The ocean model is called at evefy other 
Bate wee eeme Sisp 2O forecast prescribed 
Surface temperature for the MABL prediczlon. 
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The data sets used in this study are from the Coopéera- 
tive Fxperiment on West Coast Jceanography and Meteorclogy 
from 21 September to 12 October 1976 (CEWIIM-76) and from 5 
eom2s May 1978 (CEWCOM-78). Soe wets Corducted in the Los 
Angeies-San Diego basin. Oceanic and atnospheric data were 
collected simul*aneously throughout th2se¢ time periods. 
XBT's were deployed from tha R/V ACANIA on a regular basis 
(2-3 per day) in conjunction with atmosphere radiosondes. 
Sea surface *emperature and atmospheric data such as winds, 
temperature, 2nd relative humidity wer2 also obtained con- 
tinuously thrcugh these pericds on board the R/V ACANTA. 


Thrsé€ 24-hour periods were selected for model initiali- 


meeon and verification. Tha firs= two caseS cover a 
W8-hour period, 571970500) Oo 5721/0500 PST, ane ancy 


mueon-/8S. This péricd coincijiss with *he observazisns dis- 
@isesed in chapter 2, in which coupled air-sea interactions 


were cited as possible explanations for zthanges in both the 


flu 


atmospheric end oceanic mixes 


(D 


Leaves. Saceliic= Liegery 
showed uniformly increasing stratus coverage (thin to heavy) 


during this 48-hour period. A clear sky situation was 
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chosen for the third case. This 24—-hcur period, 10/04/1900 
mo, 1070571900 PST, was during CEWCOM-76. 


The general location of the R/V ACANIA during these 


(D 


peticds is shown in Figures 9, 10 and 11. During case 1 th 


+3 
iD 


n 


[+e 


ship was at anchor for the majority of th2 time period. 
meemcrage location is shown in Figure 9. For cases 2 and 3 


meee n/¥ ACANLA WaS Cruising slowly (2-3 kadts) returning =o 


an initial pcint approximately every 12 hours. Winds were 
primarily from the west during the cas2 1 period. The 


Winds shifted slightly to a west-northwest direction nidway 
enceugh the case 2 period. Winds for cas3 3 were primarily 
from the northwest and were much slower “han in elther case 
ior 2, 

Two predictions were analyzed for each case, one predic- 


a boundary layer noael, 


wv 


sion was based on the coupled air-s 
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and =he second vorediction was with the atmospheric no 
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With a fixed surface temperatura. Jarat2 oceanic pr 


Giction using fixed atmospheric parametsrs (persistence) was 
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that removed the wind stress and heat flux computation: 
age in the ocean model in preparation for coupling with the 


atmosphere model. 
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19-20 MAY 





meaure 9. Anchorege Position sf F/V AC 
Gtescamenaeoras Island from 5 


meeceh cas 2On OL ZeESULcS WEil be in two parts. The firs 
to.)60O6U«#ibe «presented is the coupled model's performance with 
respect tc observed oceanic and atmospheric variations. The 


second will deal with a comparison of predict2d values 


between «he coupled model and th¢ 3 


c+ 


t- 


mospharic medel. A 


{2s 


tables and figures are groupe *ogether by case number at 


me end of tne chapter. 
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Mmegure 10. Positions of R/V ACANIA from 3/18/78 to 5/25/78. 
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me nto ils FROM COUPIED MODEL 


1. CASE ji (0 
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A listing of predicted atmospheric and oceanic mixed 
layer guantities for this cas? appear in Table I. Values 
are for each 30-minute time stevo and Ssincs oceanic values 
ar2 only computed at one-hour intervals, the same vasue 
appears at two time steps. The parameters listed include: 
height of inversion (21), height ct lifting condensatior 
level (Z1lci), wind speed (Wind), atmospheric mixed layer 
potential temperature (Th), 2vaporation duct height (Zé), 
sea surface temperature (SST), ocean mixed layer depth 
PLD), and the net ocean surface heat flux (TW). 


Pred:cted values ar 


cp 


aso sown 2n F2guses “Z 
maaeough 17. BagqWec = IZ )SheCWS Ene Came evolution of predic- 
*2d sea-surface tempetatures (solid line) and the values of 
the measured surface +emperatures (dashed line). Piguce 1s 
shows the time ¢volution of the ocean nixed layer depth. 
P2gures 14-16 show ocean verticai profiles for temperature, 
moma 8 6Y Wind driven current velocities, and the salinity 
Change from an initial value of 33.5 ppt, respectively. 


Figure 17 shows predictions of 2tmospher 


\4 


= parameters. 
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High solar insulation during the daytime results in 
2 net downward heat flux (ocean warming). This can be read- 
ily observed in the net surface heat flux values (TW), in 
Table I, with the neqative sign app2aring during the daytine 
Rours. Maximum heating of th ocean occurred between 1300 


and 1400 local time. This same pattern, warming in the day- 


ct 


h 


iD 
"” 


time hours, is also cbserved in ea surface temperature 


values, Table I and Figure 12. 
The winds aprfeared to be fairly constant, 8ns"', 


between 0500 and 1100 on the nineteenth. Wind mixing coupled 


mixed 


(D 


With convective cooling at tne surface caused th 
mayer d=pth <o deepen slightly from 19.5 *9 20.9 m (Figure 
13). This deepening and net heat loss 2% the surface caused 
“he surface temperature to drop. As heating of the sea sur- 
face began tc increase due to the increase of solar radia- 
“ion, a shallow nixed layer of 14 meters formed. 

Mie wonds tap-diy decrsased from 3 0 5 ms"' at 1200. 
This coincidec with a near maximum value of tadiative hea-- 


mg calus:ng a shallow 7 meter mixed layzr to form in less 


than 1 hour. AS daytime heating cortiznzued, a shallower 
mixed layer cf 3.5 meters formed. Hit seo rOCcess 12S sSfeszr sd 


ie 
( 
iy 
py 
4 
fo 
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(D 
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fu 
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memes the ‘atternoor effect! in acoustic Li 





quite common phenomencn ir many areas of the ocean. In Fig- 
ure 14, the 7 meter mixed layer is barely discernibl2 in che 
1300 vertical temperature profile. The weak gradient at the 
base of the mixed layer only required weak mixed layer ero- 
sion (increase in winds or convective cooling) *o break down 
the shallow mixed layer. During the late afternoon, the 
winds began to increase, and with the decreasing of solar 
radiation, the shallow thermocline was ersded away. The MLD 
con*«inued deepening for the remainder of the time period dus 
@om@ind-drivenrn and convective mixing. 

Since sea-surface temperature is the sole input from 
the ocean model to the atmosphere mod2l, a comparison is 
made cf predicted and actual sea-surfaces temperature values 
(Figure 12). Since the R/V ACANIA was at anchor during the 
Majority cf this time period, this is considered the best 
cas¢ to determine how well sea surface temperature changes 
due to atmospheric forcing are peing predicted. The coupisd 
model predic*ion yielded a 90.25°9C temperature increase dur- 


ing the daytime heating period with a maxinum *emperature of 


I | 
(D 
{) 
th 


el 


su) 


15.4°C at 1600 local time. A maximum Surface temper 


ioeo°C Was observed a+ 1600 giving a 0.29 srro 
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After the daytime heating period, the surface ten- 
peratures began to drop. Th2 ninimum predicted surface tem- 
Peretuse at che 4nd cf the 24-hour period was 15.2°C which 
yields a 0.2°C temperature drop. The 0.06°C temperature 
drop prior to ocean warming nak2s the magnitude of the pre- 
dicted diurnal warming and cooling to be appreximately 
equal. The actual minimum temperature at the end of the 
period was 14.96°C. Although theres are small errors in «he 
Maegmicuds of the diurnal heating and cooling affects, the 
time evolution of this phenomenon is pradicted quite well. 


Preece Ditade10Nn WaS EOt observed dusting this period, 


therefore the surface moisture flux and surface salinicy 
flux were due entirely to evaporation. The slight predicted 


Meeeeds= in Salinity trom 33.5 to 33.5038 ppt, shown in Fig- 
ure 16, 1s due to evaporation. 

Dowie an Ve Cendi aon 2eX<isted 2n che atmosphars du-- 
mage cnis prediction period. A Swerond buoyancy flax at the 
Surface increases entrainment and should cause the inversion 


height to increase uniess subsidence (associated with shal- 


5 
ct 


Meweng St the inversicn heighs) SS Sarons wenoughmeo. Dreve 

this increase. In this case, subsidence was relatively weak 
ai 1 - : e i e e 

(-.0028 ms ) and entrainment was able to increase <=h2= inver- 


Seon height. 





Cloud occurance/dissipation was of interest because 


x 


ef its effect on radiation at the surface. mene che la ftang 
condensation level, LCL, lower than the inversicn throughout 
this time period, stratus clouds should be present during 
the entire period. This feature was observed. During this 
period the LCL decreased from 243 to 125 mneters while the 
inversion height increased from 256 to 348 meters. This 
moma =ndicace a thickening sf the clouds over this 24-hour 
period which was alsc observ2d to be trus. The air-sea 
coupling effect on predicting these cloud featur2s cannot 
be identified until the results with a constant séa surfaces 
“temperature fercing are examined. However, on the basis of 
the predicted temperature change, @ significant influence by 
ememcoupling is not expected. 
2. CASE 2 (9500 20 May £9 9500 21 May) 

A listing of predictésd atmospheric and oceanic mixed 
[ier Giantici.es for this period are given in Table III and 
syown Graphically in Figures 19 through 24. As in case 1, 
daytime heating at the surface is apparent by the negative 
meme in front of the ne= surface heat flux values. The 
Maximum downward heat flux value was not as high as in case 


Meme was probably due to thickar cloud cscsonditions for this 


ay 





case. However, increasing Surface temperatures were 


observed during «his heating period. 


4 
a 


The winds increased stsadily from 6 to 11 ms 
mieougnout the first half of thas period and slowly 


‘at the end of the period. Ths 


decreased to about 9.4 ms 
higher winds should be associated with increased wind mixin 
ct the ocean mixed layer. The increasej] winds caused ths 


mixed layer depth to deepen at the end of “this veiod. his 


also reduced the shallowing of the mixed layer during «he 


afternoon radiative heating pariod. The shallowest mixed 
layer formed at 16 meters versus 3.5 mazters in the previous 
case. In Figure 21, the 20 meter thermocline is seen to be 


very weak, and with the higher winds, th2 shallower after- 
noon mixed layer lasted only a coupl2 of hours. 

The effects of the diurnal heating and cooling peri-~ 
ods caused a 0.19°C temperature increase and a 0.15°9C temper- 
ature decrease, respectively. Both values ar2 aiso smalisr 
than those predicted in case 1. Reasons for the smallec 


predicted increase are: 1) thicker ciouds reduced <«h2 amount 


ja) 
aed 


ior oe 


a 
() 
fs 


Seecseler heating during the day, therby reducing = 


temperature increase; 2) the higher winds increased th 


ip 


Ve 


ty 


Cal mixing, causing cooler water below to be mixed into 
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*he upper layer and also increased the latent heat loss due 
29 GqaspReacdicne thus reducing the total n2t heat flux at the 
surface; and, 3) the larger alr-sea temperature difference 
increased the sensible heat loss, which also reduced the nat 
heat flux. The decrease in the downward solar radiative 
heating and increase in séensibi2 ard latent heat losses were 
somewhat balanced by the decrease in the long wave h2at loss 
due to the presence of the clouds. This decrease in the 
long wave heat loss was also the prinary reasen for the 
slightly smaller value in the nighttime cooling period. 
Comparisons cf the predicted surface temperature 
values with the actual values (Figure 19), show a noticeable 
difference. Vimo ssa baby aye te che fact that at 0700 
meceal ~2me or 20 May, the R/V ACANIA began cruising slowly 
meme tne northwest. The regional horizontal ‘temperatures 
Seructure must be taken into account. Ths southern Califor- 
mia coastai bacthymetr and associated currents greatly 


influence the sea surface tempsrature stracture. The: Csi 22= 


tt 


momma Current brings in colder water from the north as it 


qt 


heads south along ‘the coast. i@emmeor nt  CONcep=ion, the 


cr 


mibeeent begins to turn to «he south 


= 


est avay from the ccast. 


) 


BiiemeO Coastal bathymetry, Point Conception also orevides 
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sheltering for the intrusion of warm water aiong the coast 


Ss and 


@ 


from the south by the Davidson Current. The edi 
meanders of cold and warm water make the sea surface temper- 
ature structure complex in the CEWCOM experiment areas. The 
surface temperature variation observed by the R/V ACANTIA 
appears to be due to oceanographic features. 

The 19°C temperature increase in the first two hours, 
when the ship began moving t9 the northwest, indicates that 


*he ship's <‘ransit was through one of the small warm eddy 


(D 


Features. Two reasons for this are 1) Radiative heating of 
the ocean did not even begin until aftar the observed sur- 
race temperature had already b2gun to rapidly decrease, and 


2) the higher winds would tend to mix the surface waters 


with the cooler water from below, thus causing a decrease in 


a 
O 
He | 
th 


memperalure not an increase. Dre shore <ine. durat O 
this warm reéature 21so0 suggests a *transi*® “hrough an ano- 
maly. Ween whe Ship soniy Cruising at 2-3 knots, the warn 
geeMaly Stiil only lasted less than 4 hours. 


The ship turned on 2 northerly course at 1300 and 


the rapid temperature decre2ss stopped. Tt is reasonably 


ture trend (rapid decrease to slight increase) that the ship 
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mominttrally Crossing the tempereture gradient contours and 
then foilowed a parallel cours= with thes2 contcurs. This 
would make oceanic conditions more compatible with those in 
case 1. Since the 0.4°C t¢mperature increase during the 
daytime heating period compar2as favorably with the value 
observed while the ship was at anchor, the increase is 
therefore presumably due to th2 diurnal heating. 

The ship turned back to a west-northwest course 
shortly after 1800 and the surface t¢mp2rature immediately 
pegan to déecreass 1.5°C in 4 hours. Two hours after the end 
ef this pericd, 0700, the ship cane about and parallsiled ics 
previous course but in the opposite direction. An 2xamina- 


tion of the observed sea surface temperatures shows a rapid 
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=ncrease comparable to the decrease whil 
epposit= directicn. This further substantiates the observed 
large temperature changes ar¢ not due t9 atmospheric forc- 
“ng, but rather to the ship's movement through an oceeno- 
mao c eddy or front. 


Umseabl= 2t@MOSDNe=ic condizczens existed during this 


(1) 


Pesicd 3s well. The 1.5°9C drop in the air temperatures and a 
much smaller drop in the ocean temperatura (0.2°9°C) coupled 


with higher wind speeds resultad in a large surface buoyancy 
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iru x This resulted in amuch higher predicted inversion 
height change due to entrainnent. The slightly smaller sub- 
sidence rate was also a factor in the greater predicted 
increase of the inversion height. 

The lifting condensation level fluctuates up and 
down unlike tre steady decreas? in height in case 1. HOW- 
ever, the height of the LCL never incr2ased above the ini- 
tial height and so with the large increase in height of the 
mmvertsion, thickening of the stratus cloud deck is predicted 


£O continue during this period 2s wall. 


3. CASE 3 (1200 04 OCT to 1900 05 OCT) 
Redasetnag or Predicted atmespherzs and cceanic mixed 


layer quantizies are given in Table V and shown graphically 


MameeeduTes 26 thz>ough 31. Although the time period for «this 
Meme eS 2 October, there 4ar2 stili mary £satures occuriag 


mies ars Simiiar to the first two cases in May. The primary 
G@2tterences in this case, Other than being a cloud-free 
case, are 1) the air and surfaces temperature are much warmer 
end 2) a much higher value of subsidenc2? was observei. 
Daytime heating is a@pparert in predictions for the 


ear Sky Case, 2S expected. The magnituda of the net ocean 


bv 
t— 
cs 
(D 
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Seerace heat <lux is, however, smaller than the v 
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This is presum- 


the thin stratus deck condition in case 1. 
Seasonal heating rate and will 


in the 


Mocalecoupling effects. 


ably dus to the difference in 
then 


much weaker 


be considered part of the 
3 3£e 
high downward 


no * 
The winds in case 
first <=wo cases. With the light winds and 
flux during the day, a much shallower mixed 
The ocean mix2i layer depth (MLD) 
show 


surface heat 
layer might be expected. 
(also shown graphically in Figure 25) 
meter suddenly forms in less 


values in Table V 
1 

ac 0900. This shallow 

Tic =a 2c, 


began 


that an ocean mixed layer of 
than in either case 
AS the 


hour after heating 


han an 
thermocline is much stronger 


most of 


day (Figure 28). 
the shallow 


“hea 


an the afternoon, 


persists throughout 
late 

Wenas tO Mix thes 

wes Ca se. 


winds slowly increased 
deepens. 


mixed layer Without 
urface heating is much greater in 


heqech, cur 
0.6°C surface temperature increase was predicted to occur. 
Ship's movement pléeyed 


the 


As in 
Wo ciemenea On) OL 
daytime 


GOLe in 
values. This prevented the 
surface temperature increase during hs 

Sericd. Figure 11 shows the tomplex horizonc 
CEWCOM-76 experiment area for «hs 


Sesuccure in the 
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case 2, 
the large variations in observed surface temperaturs 
PESdlez ed 


a major 


Nea=22g 


al temperature 
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week of October. Peeimenis sig Ure, Peo ona sa €£O = xXxplatn 
the large 2.8°C temperature drop during the first half of 


the period as the ship appeared to be heading into warmer 


}-- 


water. The rapiG temperatur= iacrease, starting at 0900, is 


obviously due to the ship transiting into the warm eddy fea- 
ture north of Sarta Catalina Island (Figure 11). 

Unstable atmospheric conditions are also present 
mrse2=ng this time period. The high subsidence caused the 
inversion to be forced downward *o within 72 meters cf the 


surface. The lifting condensation level ditcreases evan more 


rapidly end finally becomes lower than the inversion at the 


tJ 


fast hour of the forecast period. his implies «he forma- 
om OL a iow stratus cloud deck. BeGleeLGUre 35, Log end 


stratus clouds were observed => form in the CEWCOM eéxperi- 


ment area. 


B. COMPARISON OF COUPLED AND ATMOSPHERIC MODEL 
The seccnd and final phases of this study is to compare 
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SS 
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the previous coupled air-sea boundary layer modéel cr 
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with the atmospheric model resuits using 3 constan= ss 
face temperature. Ties ee COnp2 Essen is <> cetstmane aif any 


differences existed between the atmospheric predictions, and 


j-+ 


=e SO, determine if these differances are due =90 coupled 


=o 





interactions. A listing of the atmosph=ric model predic- 
tions of atmcsphere mixed layer parameters for cases 1, 2 
and 3 are found in Tables II, IV and VI, resvectivelyv. 
Graphical presentations of thase results are shown in Fig- 
ures 18, Zo and “32. All tables and figures are grouped 
together by case number at the end of this chapter. 

The initial comparison reveals no significant difference 
<n the two predicted sets of nixed layer parameters, indi- 
cating a negliigibie effect of coupled interactions on atmncs- 
Bmeade predictions for this short time period. The 
predicted inversion heights are e¢essenti2lly identical in 
both sets of predictions for all three cases (maximum dif- 
ference of 2 meters). Mies os Sareo Lound =o be crue for 


mierda =Y and even the evaporation duct height predictions. 
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Pres ODSETVat_on 1S somewhat surprising sinc 
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nificant interactions between the su e fluxes (evapera- 


1 t) 


tion) and these parameters. The close proximity of these 
parameters ~o tne ocean surfac2 would also imply a more sig- 


Ficant impac* even with small changes in sea-surface ten- 


‘2 
j te 
4 


peratures. 
The only parameter that sesms tc be clearly affected by 
MaemecOuDling is the lifting condensation level. However, 
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even for the cloudy cenditions »f cases 1 and 2, «his shows 
only a slight deviation. A maximum difference cf only 5 
meters is observed between the predictions for these cases. 
In the clear sky case under warmer and more humid condi- 
tions, there is a 30 meter difference observed in the two 
men predicticns. However this occurred as the atmosphere 


the 


(D 


approach near neutral conditions. tes Nay  anaesat 
model's inability +o make raliabla forscasts under stable 
memospheric conditions. 

Since a comparison of predicted ocean parameters between 
a coupled and independent ocean model was neo* made, ne con- 
clusions on <his topic can be made at this time. Albenough 


tentative results indicate no definits 2dvantages using a 


iD 


Vee S 


a) 


Sempiecd model for the atmospheric predictions, the cr 


Gemeevery well be crue for ccéanic predictions. 
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On the basis of a limited study of several 24-hour pr2- 
MucteTOns , surface temperatur2? variations are not ilarge 
enough to require ocean mixed layer processes in short «ern 
atmospheric fcrecasts. Carcfully selected data sets and 
considerable interpretive efforts are essential if general 
conclusions are to be drawn on the relative roles of dynamic 
precesses in both bcundary layers. Further studies are 
required to determine the model's ability to better handle 
forecasts at longer time scales. Further studies should 
also include sensitivity analyses to define situations for 
which air-sea interaction effects are greatest. Scales 162 
PGeamlec SENSL tivity to atmospn2eric forcing are required to 
determine the impact it will nave cn acoustical forecasts. 

The CEWCOM-76 and 78 data, although able to supply 
Simultaneous and oceanic data, suffered in two major 
Tespects First, the complexity of the horizontal temperature 
structure with meanders of warn and cold eddv features had 
definite impacts on observed dsta. Second; date was col= 


iscted while the R/V ACANIA was transicing through a ccmplex 


G2 





Seearmecetiperac-Ure Structure. Changes in observed conditions 
under these circumstances can and did mask the interaction 
processes of interest. 

Although the primary Best ie Gf this study indicate 
coupled air-sea models are probably not necessary f9r short 
“erm atmospheric forecasts, a coupled mod2l still has advan- 
tages, e¢specially if long tern forecasts are to become pos- 
sible. Such a modél would not be the single one-dimensional 
one considered here. However, this model offers many pcoeten- 
tial benefits as a multi-purpose model wit Single station 
assessmert and prediction capabilities in the absence of 
shore-based support. A scund velocity profile (SVP) predic- 
meen Ccapabil=* GOuuld ~~ bemeasi dy ~incorpotated 2nto this 


mem.  ii2n Conjunction with prasent predistion capabilities, 


4 multi-purpose nodel would bé available For a wide range of 
needs. wiccowmm=xed leyersdepch “and SVP predictions car be 
Meem in predicting ASW sensor capabilitizs. Environmental 
Ssrfeczs on tne perccrmance of weapons and sensors that 


cepend on eiectreomagnetic and 2ilectro-optical wave propada- 


=icn can be evaluated using <h2 same nodel. 
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Figure 12. Time Evolution of the Sea-Surface Temperaturs-. 
Case 1. Predict2d values are shown oy the solid 
+7ne. Opserved values are shown by the dashea 
Pale. 
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Figure 13. Time eee eo” of Jcoean Mixed Layer Depth (MLD)- 
Case 1. 
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Figure 14. Selecteé Ocean Vartical Temperatures Profiles- 
Case 1. 
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Figure 15. Selected Ocean Vertical Profiles of the 
DO==zZonec wend DE=yer Curcsces- Case 1. Sast- 
west compcnent u (solid lines) NGET Rn =SCuch 


component v (dashed). 
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Figure 18. Graphical Display of the MABL, Model Predicted 
Atmospheric Mixed Layer Quantities- Cease 1. 
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Figure 20. Time Evolution of 3cean Mixed La 
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meggre 26. Time Evolution of the Ssa-Surfiacs Temperature-.._ 
Case 3. Fredicted values are shown by the soiid 
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Figure 27. bone Evolution of Jcean Mixed Layer Depth (MLD) - 
ase 3. 
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Figure 28. Selected Ocean Vertical Temperature Profiles- 
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Figure 29. Selected Ocean Vertical Profiles 
Horizontal Wind Driven Currents- 
west compcnent u (solid line) No 
component v (dashed). 


31 





Proftle times 


338 9e8 1438 
8 AT CUUT PUTT P TTY TTP TOT TET TNT | PPTL yTT 
& = ‘ 
» 
o 
& -22a 
-28 . 
—25 8 25 58 75-25 8 25 58 75-25 8 25 58 75 
Salem cyveecmancge: X 1O00H Coot) 
Pegure 30. Selected Change in Jceai Vertical Sealine ey 


Profiles- Casé 3. 


82 





qn (g/kg) 








Height Ckm) 


38 


22 


Cm) 


me 18 


Pesanize 31. eee yee 9f +he Couoled Model Pr 
Atmospheric Mixed Lay 


SO) ee 8. 


Time iS 





Time 13 


Ges = G 
@ 


= 
Sr Quantitzies~- Case 3 


83 





18 -3 





er 
iv 
N Pe eee 
Cee 
= 
Ss 
845 elegans 


Time 


o~ 
£ 
a 
ww 
» 
Soy 
om 
0 
“as 





38 
S 23 
Sees? 
ar 18 
8 pee et la ye 
‘9 Time . 


Peme@eieeee Ss 2, Graphical Display of “he MABL. Model Predicted 
Atmospheric Mixed Layer Quantitiss- Case 3. 


84 





*potisgd uctzyoOtpsid 
¢ asey butang suotitsod pnc{ta snz%azg§ pue bog vee aanbta 


wet Ui Wert a | 


a et | ee 
$12 99519 / 


Cte abe 
c ’ 











‘2 6'U 

Orr etd ; sO Va Ser asaice 
“4, be | 

ONL a OFS ae AAAS ©7000 vorgunt 






4 
6°12 §S¢ ao ae 
ct < -.on VW wy oS 
mean A z'2t))O¥GA \ 
‘ 09 Lion) ’ 








obpy ye Gos 






a ae ~ 
v% P= ~~) are 
> tb 
Ne ql) 
N 
N 
ee “* 
ed aN aa 
"HY } Be 
ae 
wo” 
: yo FD u's “tues 
(4.0 
i'ue seats 
92 Wooma0 wey y 
A 
Cd 27 Gas Gog 
(Opn Gg Aagponns; e 





Gabry pserveces ggevd.aey UIE BIO MPD ° 9D 


J TAN-SdN . oTtantuy °°) 


85 





hse lene ERENCES 


Bees, Dasa lov? se cebea Ana aA Ssociated Offshore Fog: 
Forecasting with a eases o MSCere. Mes. Lhes*s, Naval 
Peeteaoeddua-zs School, Monzerey, CA., 112 pp. 


Beower, D.A., 1992; “Marine Atnospheric Boundary Layer and 
ipversior Forecast eee M.S. Thesis, Naval Postaraduate 


penool, Monterey, CA., 136 pp. 

Deweadson, K.L., C.W. Fairall, P. J. Boyle, 2nd G.E. 
pempacner, 1983; "Versticaeti:on Of ac AtmaSs0heric Mixed-Layer 
Meme. £92 a Coascal Reg:on", submitted,CA., 39 pp. J. Appl. 
Meteor., 28 pp. 


Wo reatecionh OF Convective. Mixed La 
st 
43 


Deardorf, J.W., 1979 
Entrainment for Real 
AzMOS. Sci-, 36, 42k 


Vv 
se Gamba ng Invertsicn Stzucture", J 


~-L. Davidson, Tf 
"Atmospheric T 
N 


K 
ng CEWCOM-76", 
noo, Monterey, 


S in 


ence, ha tan Dave ison, and 


Bas rali, C Send . n 
Pics Mieese Macone Boundary L 
Vie o 
= 00d 


Peudehan, 1 
Measurement 
SaacON-/8", 
Monterey, Ca 


Pearveeem someN? Colas 1s 
aionuwui, Navat Postgraduate § 


Moreeell, C.W., K-L. Davidson, 2nd G.E. Schacher, 1981: "A 
Review 2nd rvaluation or Integrated Atmoshperic Bourndary- 
[eee cevodcis", NPS-63-81-004, Waval Postgraduate School, 
Memecersy, CA., 8Y pp. 


Cn ake ee eA Teer SD duction <o 
hysics." New York, New York: Academic 


re 


Seo cee Alc Cam Ma x 
See ic eeu a cos 


an, J-P., 1982: "Single-Station Assessments of the 
BemoGate TOrC nq On tas Maaine ,Acmospheric Boundary 
meio. Thesis, Naval Postgraduate School, Mont= rey, 


5§ DD. 
mieoy, D.M. andG.R. Seckel, 1978: “Large Scale Air-Ssa 
Mbeeeee ONS at Ocean Station 7", NOAA Tech. Rep. UMPS 
SSRF-696, 44 pp. 8 as 


Meevesel, T., 19605 "Factors Atlecting the T[emperature of 
Surface Layer of the Ocean", Soc. Scieni. Femica. Ccoumenc. 
Paysico.-Mathem., 25(1), 1-136. 


—_—— ke ee ee ee ee 


86 





(feelers mats, 197558 Ine Selinitty EEfects in a Mixed Layer 
Ocean hodei"' Umea OGScdIoage. , 6, 29-35. 


kel, G.R. and F.. Beaudry 
Sky over the North Pacifti 
phys. Union, 54, 1114. 


Sioa the Rad: action £5om Sur 
Coane, SOS,8 Lears. Al. 


a cae 


, 1 
co 


oO (mimo, occ the One—-Dimensionality of the ED RSE 
Ocean Mixing and the Role of Aivection during the 3o 
Experimnen.", MM. 3 3 Thesis, Naval Postgraduate Schocl, 
Monterey, ch., Ppp. 


Peaetan, elo vce soca SUrtAcCS Thermal Characteristics in 
the Vicinity of Southern Californias Channel sirandas' 2 
Department of Oceanography, Naval Postgraduate School, 
Pemecrey, CA., 61 pp. 


Pengo, N., oS. Nicholis, and C.L. Wrench, 1982: "A Field 
Beudy Of Nocturnel Stratocumulus: ate iaeerd 2 ‘gh Resolu-ion 
Radiative and Microphysicai Observations", Quart. J. R. Met. 
Bac. ; S, ldo—leo. 


midges owA. ahd Wel. Buscnger, 1981: “A Model for 
eMiraa son ent TReCa Soule Topped Marine Boundary Layet- Par 
Be Medel Discription and Appi isa sion orm "eo lara: One peeee 
Foisode", we A2mes. Sci. , Z2AS=2 229s 


stage, S.A. and Vote eicende:, 1981: "A Model for 
Pntrainment into a Cleud- -Toppel Marine Boundary Layer- Part 
nL: Discussior of Model Behavior and Comp2rzison with other 
Pomeis", J. Atmos. Sci., 38, 2230-2242. 


mace os, I!704s "NR Study of the Main Physical fF 
Governing the Qesanogzaphic Cond aons Of Station 
Northwest Pacific Océan Ph. D. Thesis, Jniversi 
mexkyo, 264 pp. 


Geagerss, H. end A.G.M. Dreidonks, 1981: "Basic Entrainmen 
Equations for the Atmcespheric Boundary Layer", Boundary 
Layer Meteorcigqy, 20, 515-531. 


87 





‘oO 


10. 


11. 


Mert aa DISTRIBUTION List 


Defense Technteal © 1fOrmation Center 


Cameron Station 
Aiexandria, VA 22314 


Library, Code 01742 
Naval Postgqradua*e School 
Monterey, CA 93940 


PEMreasScre NRedemaenera, Cod2 63Rd 
Naval Pcstqraduate School 
Monterey, CA 93940 


Professcer C.N.K. Mooers, Code 68Mr 
Naval Postgraduate Schoo. 
Monterey, CA 93940 


Professor K.L. Davidson, Code 63Ds 
Navel Postgraduate School 
Monterey, CA 40 


Professor R.W. Garwood, Coie 68Gdad 
Navel Postaqraduaze Schoo 
Monterey, CA 93940 


eee ne -wccende net ecas 35 3'Ga 
Naval Postgraduate Schcol 
Monterey, CA 93940 


DEsCe tor 

Nava. Oceanograp 
Naval Observator 
34TH and Massach 
Weshington, D.C. 


Ey Divison 
y 
u 


NOW 


Commander 

Naval Oceancgzrtaphy Command 
Hole o. LOG 

Boweom, LOUutS, MS 39522 


Commanding Officer 
WeveemOcecancgqraphic Offic: 
NSUL§sStation 

Bevaeoe.s LOUTS, MS 39522 


Commanding Officer 


Fliset eee aoe Center 


Non =erey, Ca 


88 


Copié¢s 


10 





Commanding Officer . . Se 
Naval Ocean Research and Development Activity 
Nol terocad LOR . 

Bay) cee Lous 5 ts 39522 


Commanding Officer we oo 
Naval Environmental Prediction Research Facility 
Monterey, CA 93940 


Chairman Oceanography Department 
U.S. Naval Acaden 
Anapolis, MD 2140 


Gime wOr Naval hese 
800 N. Quincy Stree’ 
Arla ng@ecn, VA 22217 


Office of Naval Research (Code 480) 

Naval Oceanography Research and Development 
ACTIVITY 

lott stat 1on 

Beyrot.e LOWes, DS 39522 


Setcn 1 eee gera. Scr Offics 

Office of Néval Research 

eo eS Institute of Oceanography 
La Jclla, CA 92037 


LupEary | 
Bere RS one 
Pom. SOx 236 
La Jolia, CA 


tute of Oceanography 
B70 30. 


~Jt? 


ies a Ly 

Department Of Oceanography 
peace Creer ang Or 
Seattic, WA 98105 


03 
Senmeys.aro, CA 92073 


Leary _ 

SEneet Of Oceanography 

Oregon State University 
Corvallis, OR 97331 
Commander. 
Oceanographic Systems Pacific 
BOX Q 


Pearl derbor, BRI 96860 


89 





ceo 


24, 


ZO; 


Commanding Officer 
eer: 
CDR D, Hinsman 
Dr, A, Weinstein 
Naval Environmental Prediction Research Facility 
Monterey, CA ; 93940 


CDR J.J. Jensen, USN, Code 554731 
Office of Assistant Secretary of the Navy 
Assistant for Environmental Science 
Washington, DC 20360 


Commanding Officer 


Attn: 

LT T. McPherson 
COMNAVSOUTH 
BOX 168 


FPO NEW YORK 09524 


Commanding Officer 
ttn: 
eae anor ©. Lowe inlan 
Naval Eastern Oceanography Center 
Naval Air Station 
MOrtolk, VA 2351) 


90 














200250 


Thesis 
0473 O'Loughlin 
e..) Formulation of a 


prototype coupled 
atmospheric and 
oceanic boundary lay- 
er model. 





